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Abstract 



^ \ Nonperturbative flow equations within an effective linear sigma model cou- 

pled to constituent quarks for two quark flavors are derived and solved. A 
heat kernel regularization is employed for a renormalization group improved 
effective potential. We determine the initial values of the coupling constants 
in the effective potential at zero temperature. Solving the evolution equations 
with the same initial values at finite temperature in the chiral limit, we find a 
second order phase transition at T c 150 MeV. Due to the smooth decoupling 
of massive modes, we can directly link the low-temperature four-dimensional 
theory to the three-dimensional high-temperature theory. We calculate the 
equation of state in the chiral limit and for finite pion masses and determine 
universal critical exponents. 
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1 Introduction 



At zero temperature and zero chemical potential, the chiral symmetry of Quantum 
Chromodynamics (QCD) is spontaneously broken. One expects that the chiral sym- 
metry will be restored at sufficiently high temperature and a phase transition occurs 
separating the low temperature and high temperature regions. This phenomenon 
may be realized in high energy heavy-ion collisions at RHIC (Brookhaven) and LHC 
(Geneva) and is therefore a subject of general interest. It is a matter of intense de- 
bate at what temperatures the phase transition would occur and what its nature 
is. 

In this paper we study the chiral phase transition at finite temperature both in 
the chiral limit and for finite pion masses. We present flow equations which link the 
universal critical behaviour at finite temperature to the well known physics at T = 0. 
Calculating the effective potential in the framework of the renormalization group 
theory with a heat kernel regularization, we investigate the critical behaviour of the 
chiral phase transition for two light quark flavors and calculate the critical exponents. 
It is still an open question whether full QCD with two massless quarksQ and the 
three-dimensional 0(4) Heisenberg model belong to the same universality class. We 
use the linear sigma model with quarks which at T = exhibits a spontaneous 
breaking of the 0(4) symmetry to an 0(3) symmetry. We show how the 0(4) 
symmetry is restored at sufficiently high temperature. Lattice investigations hint 
that the two flavor chiral phase transition lies in the 0(4) universality class |1] from 
which the behaviour of the condensates and the long distance properties at the 
critical temperature ensue. If this is confirmed the transition would be basically 
driven by the pions and the chiral sigma particle. 

It was realized that a perturbative expansion of the effective potential breaks 
down near the critical temperature due to infrared divergences and the nature of 
the phase transition cannot be investigated in any finite order in perturbation the- 
ory. The situation was partly resolved in by means of the average potential which 
employs renormalization group ideas. Constructing evolution equations for the de- 
pendence of the effective potential on a variable infrared cutoff this approach sums 
effectively the relevant class of perturbative diagrams in the complicated infrared 
regime. We follow this method and calculate in a nonperturbative way the effective 
potential using a variable heat kernel cutoff. The advantage of this cutoff method 
is that the evolution equations can be obtained in analytical form which makes the 
physics content very transparent. The aim of this paper is twofold: We demonstrate 
how the renormalization group improvement of the mass and coupling constant evo- 
lution give the main physics of the chiral phase transition, i. e. the equation of state 
away from the critical temperature. At low temperatures one sees the success of the 
model handling the pion and sigma degrees of freedom in a comparable accuracy 
to chiral perturbation theory. At intermediate temperatures the quarks take over 
as main agents contributing significantly to the change of the chiral condensate. 
The high temperature behavior of the linear sigma model has not enough degrees 
of freedom to model full QCD lacking in particular gluons. 

By explicitly constructing the thermodynamic equation of state we think a more 
realistic assessment of the contents of the model is possible. We also think that our 

lr The influence of the strange quark is neglected in this work. 
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method is very pedagogical making this powerful method available to physicists who 
are interested in an intuitive access to an otherwise very technical subject. 

The paper is organized as follows: In section ^ a general summary of our for- 
malism for an iV-component 4 theory coupled to Nf quarks is given. After the 
calculation of the effective potential our ansatz for renormalization group improve- 
ment is presented. In section |3] the nature of the chiral phase transition is elucidated. 
The critical behaviour is discussed and critical exponents are determined. In this 
section the free energy density i.e. the thermodynamic potential for the transition 
and the equation of state is calculated in the chiral limit. Section |] is reserved 
for a discussion of the same properties for finite pion mass. Section [| presents our 
conclusions. In order to improve the readability of this paper we postpone extensive 
equations to appendices. 



2 Flow Equations 

Wilson's original renormalization group ideas have been used widely in perturba- 
tive and nonperturbative calculations. A nonperturbative approach has been used 
on lattices and also implemented in continuum Quantum Field Theory by Wegner, 
Houghton and Polchinski 0,(0 • Flow equations describe the average of an effective 
action over a volume and are the continuum analogue of a block spin transformation 
on the lattice. All degrees of freedom with momenta larger than an infrared scale 
k, the coarse graining scale, are effectively integrated out. The flow equations are 
ultraviolet and infrared finite through the introduction of an infrared cutoff function 
fk, which obeys several requirements described in the next subsection. The solution 
of these equations provides the full effective potential, once the infrared cutoff k is 
removed (k — > 0). On the other side in the limit k —>■ oo or k equal to some ultravi- 
olet cutoff A the effective average potential results in the classical potential where 
no fluctuations are included. This ultraviolet cutoff may be associated with the 
highest momentum scale for which the theory parameterizes the physics adequately. 
All quantum fluctuations with momenta larger than A are included in the poten- 
tial at this scale, thus no further necessity for an ultraviolet regularization scheme 
arises. The effective couplings at the large scale A specify the potential completely. 
The knowledge of the /c-dependence of the effective average potential then allows 
to interpolate from the classical 'bare' potential at k = A down to the effective 
'renormalized' potential for k —>■ including more and more quantum fluctuations 
as the evolution proceeds towards zero. 

This method of the flow equation has been developed intensively by Wetterich 
and his collaborators [0. One important ingredient of the evolution equation ap- 
proach by this group is the appearance of the exact inverse propagator in the flow 
equation, which makes the RG-improved one-loop equation exact and intrinsically 
complicated ||21|| . The solved truncation of this exact equation is very similar to the 
RG-improved one loop equations used here where the bare couplings are replaced by 
their corresponding full fc-dependent expressions. The technically difficult part in 
the exact RG-approach is not the establishment of an exact relation, but rather the 
choice of a suitable nonperturbative truncation scheme in order to solve this highly 
nonlinear coupled set of differential equations. We use a simple truncation scheme 
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and analytical flow equations can be derived. We believe that a deeper connection 
between Wetterich's and our approach should be possible. 

The flow equations are nonperturbative, i.e. also applicable for problems with 
large couplings, loosely spoken, because the evolution of the different couplings is 
interlinked and no particular truncation in coupling constant is used. The solutions 
of the flow equations give an accurate picture of the high temperature phase tran- 
sition in four dimensions. It is possible to study the phase transition at the critical 
temperature and to calculate the critical exponents with surprisingly high accuracy. 

2.1 The Effective Potential 

In this subsection we develop the formalism starting with the derivation of the flow 
equations for the linear SU(2) x SU(2) sigma model combined with two quark 
flavours. The linear sigma model is a hybrid model which contains quark and 
mesonic a- and 7r-fields as collective degrees of freedom in order to simplify the 
strong qq interaction. The cutoff A plays the role of a compositeness scale, below 
which the effective o- and 7r-fields can be treated as elementary. For the moment 
we leave the connection of this theory to QCD as an open question and see how far 
we get with this model. Preliminary investigations in this direction have been made 
in ref. M. 

The partition function at zero temperature is given by 

Z[J = Q] = J VqVqVaVn exp{- J d 4 x(C F + C B )} . (1) 

Here we omit any external sources J, because in the one loop approximation the 
effective potential and the generating functional for one particle irreducible Green 
functions do not differ. The Euclidean^] Lagrangian in d = 4 dimensions looks like 

C F = q(x) {^ E d E + 9 + «T7f7 5 )) q(x) , (2) 
C B = \((d,af + (d^) + ^(a> + ^) + ^ + nr-c* . (3) 

The parameter c is an explicit symmetry breaking term in cr-direction, which gives 
the Goldstone boson a finite mass. In order to avoid later any recurrences we derive 
the effective potential with an explicit symmetry breaking term c. The chiral limit 
can be extracted by putting c to zero in all equations. 

Taking advantage of the properties of the minimum O of the potential we can 
eliminate the mass parameter ml in favour of Ao and <pQ. Then the potential can be 
written with m\ = c/(po — ^o<t>o as 

V^-tf)'-^-. („-!-) (4) 

where we introduce the 0(4) vector <fi = (a, tx) as a shorthand notation for the meson 
fields. 

2 We denote all Euclidean operators with an index 'E'. 
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Formal integration over the fermions yields a non-local determinant, which can 
be regularized by the heat kernel representation. We use the abbreviations M(x) = 
a(x) + iT7f(x)-f 5 , (MM + = (f) 2 ) and D = ^ E d E + gM(x). 

Z[J = 0] = J PaPvrdet (-? E d E + gM(x))exp{- J d 4 xC B } 

= J VaVnexp^Tr log DD + - J d 4 xC B } (5) 

= J VaVn exp{— / d 4 xC B + 

-- f — [d 4 xtr(x\e-<- d e + ° 2MM++ ^ dM+ ^\ x)} . 

We remark that the lower boundary of the integral over the proper time r reflects the 
ultraviolet cutoff scale A, while the upper limit of the proper time integral represents 
the infrared scale. Later, we will choose an infrared cutoff function in such a way 
that no ultraviolet cutoff in the proper time integration is needed. Details concerning 



the heat kernel representation and definitions can be found in ref. |T2| . 

Because we are interested in the calculation of the effective potential no wave 
function renormalization corrections are taken into account. This means that we 
set the wave function renormalization constant Z k — 1, which simplifies further 
calculations because derivatives resulting from the heat kernel expression can be 
omitted. This is not a principal restriction of the method. On the contrary the 
chosen heat kernel regularization is designed for a derivative expansion and higher 
derivatives can be taken into account not only for flat manifolds but also for curved 
manifolds 

We find for the partition function at zero temperature the following expression 
which is exact up to the omission of derivatives <9 M in the fermion exponent 



Z\J = 0] = / VaVnexp{ — f°° — [ d 4 xtr(x\e-<- d ^ +92 ^\x} - [ d 4 xC B } . 
J 2 Ji/a 2 r J J 

(6) 

The fermion determinant is in general part of the functional integration over the 
remaining meson fields weighted by the meson Lagrangian C B . In the following we 
limit this integration over the fermions to the one-loop level. Thus we decouple 
the effect of modified fermions back on the meson dynamics, pull the non-local 
determinant in front of the meson integration and replace the meson fields <p 2 with 
the vacuum expectation value (VEV) 4>o- The total effective action splits into two 
pieces and becomes a sum of fermionic and bosonic terms each depending on </> : 

- \nz = r(0 o ) = r F (0 o ) + r B (0 o ). (7) 

In a plane wave basis for the diagonal part of the heat kernel \T2\ the effective action 
for the fermions T F has the form 

where the trace runs over color-, flavor- and spin-space and gives a factor 8iV c for 
two quark flavours. 
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The remaining meson integration is also performed on a one loop level and is 
obtained by a saddle point approximation. The classical mesonic potential Vq = 
mQ0 2 /2 + Ao(0 2 ) 2 /4 — co is expanded around the minimum, which corresponds to 
the VEV 0o- Thus taking quadratic fluctuations of the mesons into account and 
using again a heat kernel regularization for the fluctuation determinant one finds 
the effective action for the bosons: 

,b,, , i t m. t°° dT r A L "(» 2 +iiS£->' 



r«(0„) = --y^y^-y^i^e— v| . (9) 

Here the trace runs over the (4 x 4) fluctuation matrices in (a, 7?)-space. Note the 
opposite sign of the bosonic and fermionic effective actions. In the cx-model the 
calculation of the second derivatives of the potential Vq is straightforward and the 
trace can be evaluated explicitly for N(= 4) components. One has to determine 
the eigenvalues of the fluctuation matrix evaluated at the minimum of the potential 
i. e. the masses of the system. With an explicit symmetry breaking term c the trace 
is given by 



tr N e 



The first term on the r.h.s. of eq. (|I0| ) represents the mass squared eigenval- 
ues associated with the three pions while the other term describes the cx-meson 
(cf. eqs. (fj9|)-(|6lD). In general this trace is a sum over all exponentiated squares of 
mass eigenvalues times the proper time r. 



2.2 The Renormalization Group Improvement 

The effective action r(0 o ), eqs. (0)-(0), contains so far only fluctuations up to the 
one loop level and depends on the ultraviolet heat kernel cutoff A via the proper 
time integration. In order to cut off the fluctuations in the infrared we introduce 
into the proper time integrand a universal fc-dependent function 

fk(r) = g(x) (11) 
x = rk 2 (12) 

with a dimensionless argument x = rk 2 which depends on the product of the proper 
time and the square of the infrared cutoff k. This function has to satisfy three 
conditions: 

We require that g(x — > 0) = fkijk 2 — > 0) = 1, since the action Tk with infrared 
cutoff should tend to the full effective action r at k = 0, i. e. in the limit k — > the 
infrared cutoff is removed. 

The function g must satisfy g(x — > oo) — > 0, to suppress modes with fixed "virtu- 
ality" k for large r. Remember that the upper limit of the proper time integration 
represents the infrared regime of the effective action which is regularized by this 
condition of the function g. 

Thirdly, the first derivative of the cutoff function g has to obey 

g \ x ) = - x 2 h(x) (13) 
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with h(x) being any regular function in the vicinity of the origin. This condition 
assures that the renormalization group equations for the /c-dependent effective action 
are ultraviolet finite. The cutoff function regularizes the derivative of the effective 
potential with respect to the infrared parameter k, therefore a further ultraviolet 
heat kernel cutoff A is unwanted. Due the above condition the ultraviolet divergences 
in the effective action itself are not affected by the presence of the infrared cutoff 



21] , |2o[ . It is this condition which distinguishes the cutoff function from a simple 
minded mass cutoff. 

One possible choice for the function g(x) which satisfy all these requirements is 

g (x) = e~ x (l +x + ^x 2 ) (14) 

The function g depicts a smooth cutoff of the momentum modes and allows for an 
elegant and simple form of the evolution equations. 

A smooth cutoff is unavoidable if one wants to have a well defined derivative 
expansion. For the effective potential alone also a sharp cutoff would do 0,0.0. 

Even in this case the above choice has great advantages over a sharp cutoff, at 
least in the proper-time integral formalism where it may cause unpleasant oscillations 



(cf. e.g. [28 



The effective action with the cutoff function /fe(r) has now the following form: 

r k [<f) ] = fd A xV k (<p ) (is) 



with 

V-^-i/fAM/jgjj^e-^^'-tre-^^)} . (16) 

This equation defines the one loop effective potential with an additional infrared 
cutoff function. In a further approximation we truncate the effective potential at 
terms oc (0 2 ) 2 . Such a procedure leads to simple evolution differential equations for 
few couplings (A, m 2 ). These evolution equations are improved by replacing the bare 
couplings (A , ml) and the constant VEV (0 O ) with the running couplings (A^m 2 ,) 
and the running VEV (4>k)- This renormalization group improvement resums the 
relevant infrared divergent Feynman diagrams and gives the important nonpertur- 
bative evolution of the potential. The substitution can be compared with ordinary 
renormalization of n-point functions in perturbation theory, where one replaces on 
the one loop level the bare couplings with the renormalized ones. This simple strat- 
egy is pursued in this paper and has many advantages. Another possibility is to 
derive an evolution equation for the full potential using the heat kernel regulariza- 
tion. This is postponed to another paper in the future |16 . 



The evolution of the potential itself with virtuality k leads to two phases and 
therefore two sets of equations. In the chiral limit the minimum of the potential 
(VEV) will evolve from a large value to zero in a second order phase transition. For 
finite pion masses (c 7^ 0) this transition becomes a crossover and the two sets of 
flow equations merge to one set of equations, as we will show analytically. In order 
to simplify and clarify the following discussion we start with the chiral limit and set 
c = for the rest of this section. 
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For high virtuality k the system is in the symmetric phase which is defined by 
a vanishing VEV (<f) k = 0) of the potential. Through the derivative with respect to 
the fields 

* == w (17) 

we can define the mass ml and the coupling constant A& in the symmetric regime 



mi 



= 4>l = o) 



y := W* = $ = 0) • (18) 

The effective potential at the minimum is given by 

v k :=V k $ = <j> 2 k = 0) . (19) 

Taking the derivatives of equation (103) with respect to the scale k we find the 
following coupled sets of flow equations for the symmetric phase (0| = 0): 

kd\ k _ avm 
2 dk ~ k dk ■ 

In the spontaneously broken region the VEV is finite {<f> k ^ 0) and the mass 
parameter ml tends to negative values. We prefer to parameterize the evolution of 
the potential in this region in terms of A^ and the minimum of the potential <p k , 
which is defined by 

YLfa) = . (23) 

This equation enables us to find the evolution of the minimum <p k . By again taking 
the derivatives with respect to the scale k we find the flow equation for <p\ in the 
broken phase (4>\ 7^ 0) 



dv k , dV k 



k^ = k^ , (24) 

kd£_ k dvi^l) k avKrt) /o-x 

2 dk 2V£'($) dk X k dk { ' 

kdX, _ W 

2 dk - k 8k ' lJ ^ - U • (2b) 

Here an important remark must be added, that we first evaluate the derivatives 
on the right hand sides for constant couplings A , m\ and constant VEV (f> with 
respect to k. Thus the scale derivative acts only on the infrared cutoff function f k {j) 
inside the proper time integrand. 
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3 Equation of state in the chiral limit 



In this section we discuss the flow equations for the two flavour linear sigma model 
in the chiral limit. We present the numerical solution of the flow equations with 
respect to the scale k and generalize the equations to finite temperature. 

3.1 The evolution for T = 

In the symmetric phase at zero temperature we use the equations (|20|j2l|j22]) and 
get the evolution equations for the free energy density, mass parameter and quartic 
coupling. Due to the choice of the cutoff function /^(r) the resulting equations are 
very transparent. The heat kernel cutoff is not a simple mass cutoff, in effect the 
derivative cf. eq. (14) of the cutoff function is simpler than the cutoff function itself, 



therefore the evolution equations have such a concise form: 



dvk _ 20 

k dk ~ 2(4vr)2 fc ' {I n 

kdm\ 3X k k 2 1 AN c g 2 k 2 

H / . NO 5 



2 dk ~ (4tt) 2 (1 + m 2 k /k 2 ) 2 (4tt) 2 

(29) 



kd\ k 12X 2 k 1 8N c9 4 



2 dk (4tt) 2 (1 + ml/k 2 Y (4tt) : 



The factor 20 in eq. (|27| ) arises from the combination 4 — 8N C (4 mesons minus 
4NfN c = 8N C fermionic degrees of freedom). In the symmetric phase the pions are 
degenerate with the a-meson and the squared mass parameter is always positive. 
The VEV <pk vanishes and does not appear on the right hand sides of the evolution 
equations. In the broken phase the pion mass vanishes in the chiral limit. The 
constituent quark mass is proportional to g<j) k and is zero in the symmetric phase. 

To solve the evolution equations we need three initial values of v k , X k and m| 
at k = A in the UV. We start the evolution in the symmetric phase deep in the 
ultraviolet region of QCD at A = 1.2 GeV and set v k = 0. The other two parameters 
Xk and ml are adjusted to get (fi k =o = 93MeV = f w at k = and the transition to 
chiral symmetry breaking at k — k x sB = 0.8 GeV. This leads to m 2 k=k = (550MeV) 2 
and Afc=A = 40. 

In addition we need the Yukawa coupling g which we fix to g = 3.2 for all k 
in order to get a constituent quark mass of 300 MeV at the end of the evolution. 
Because we do not take any wave function renormalization into account g has only a 
weak scale dependence. In the work of Wetterich et al. [§] the main contribution to 
the strong evolution of the Yukawa coupling stems from the evolution of the meson 
wave function renormalization. 

In fact, the initial value for the quartic coupling Afc=A turns out to be not very 
important due to the Gaussian fix point X(k — > 0) — > 0. This behaviour can be seen 
in figure [I], where the evolution of the quartic coupling towards zero for different 
initial values is plotted. Near k = A the evolution of X k is dominated by the 



bosonic part, the first term on the r.h.s. of eq. (|29D whereas the evolution of m% 



responds to the fermion loops. When we reach chiral symmetry breaking at the scale 
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1200 



k [MeV] 

Figure 1: The evolution of with respect to the scale k for 
different initial values. 



k = k x sB ~ 0.8 GeV, we switch to the equations for the broken phase characterized 
by + 0: 



k 



dk 
kdtf, 
2 dk 

k d\k 
2~dk 



k 4 
2(4tt) 2 

3k 2 
2(4tt) 2 

(4tt) 2 



3 + 
1 + 
1 + 



1 

;i + 2A.02/P) 
1 

3 

(l + 2A fc 0l/A; 2 ) 3 J " Jm 



8N C 
(1 + g^l/k* 
AN C k 2 g 2 
~ (4tt) 2 X k 



1 



iiV, 



V 



1 



(30) 
(31) 
(32) 



The first two terms on the r.h.s. of the above equations are again related to the 
boson fluctuations while the last parts are connected to the fermionic contributions 
to the evolution. 

The constant factor in the first part of the upper equations comes from the three 
massless pions while the second part of the bosonic term contains the contribution of 
the cr-meson with squared mass In the broken phase the cr-meson is massive 

(4>k 7^ 0) and the a-meson contribution decouples from the evolution for small k 
values in the IR region, when k <^ 2\ k (pl; whereas the massless pions contribute 
permanently to the flow in the chiral limit. The fermionic quark part of the flow 
equation enters with the opposite sign to the evolution and the quarks become 
more and more massive in the infrared. The factors have a similar origin as in the 
symmetric phase. The contribution from the pions is three times larger than the 
cr-meson contribution. Two flavors , two particles and anti-particles and two spin 
degrees of freedom yield a factor of 8iV c . 
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Figure 2: The evolution of the minimum of the effective potential 
(fi k for k < k x sB and the evolution of the mass parameter for 
k > k xSB . 



The functions in squared brackets can be called threshold functions, because 
they describe the smooth decoupling of massive modes from the evolution towards 
the infrared limit of the theory given by the scale k. With the heat kernel method 
and our choice of the function /^(r) we can calculate these threshold functions 
analytically. They are of the form (1 + mass 2 //c 2 ) _p where p is some positive integer 
power. So the quark loop gives a factor (1 + g 2 <t>l/k 2 )~ p with positive p which 
switches off the quark contribution to the evolution of the couplings for k <C g<frk- 
Similarly, the square of the sigma mass over the infrared scale /k 2 suppresses 

the massive boson loop. At small evolution parameters only zero mass pion loops 
drive the evolution. This is very intuitive and coincides with the ideas of chiral 
perturbation theory. In fact, it has been also shown p0| that this model reproduces 
the chiral expansions even beyond leading order. As we will see in the equation of 
state the suppression of the quark terms for the evolution does not mean that the 
quarks do not contribute to the thermodynamic functions like the pressure or energy 
density below the phase transition temperature. In this respect the model has the 
same defects as more naive NJL models without quark confinement. 

We remark that in the heat kernel expression for the effective potential the inverse 
fermion propagator enters quadratically in the combination DD + therefore we can 
use the same infrared cutoff function as in the bosonic integral without breaking 
chiral symmetry. If we want to take into account the running of the Yukawa coupling 
g the fermion propagators also have to be squared. One also recognizes the signs of 
the bosonic and fermionic contributions to the /3-function. The bosons lead to an 
infrared stable (ultraviolet unstable) A-coupling, whereas the fermions counteract 
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this tendency. Going from high k to low k one sees that the mesonic self-interaction 
Afe will balance at intermediate values of k, whereas in the far infrared the boson 
term wins (cf. figure ||). The chiral symmetry changes at k = k x sB ~ 800 MeV, 
where the system goes over from the symmetric phase to the spontaneously broken 
phase. At this scale k eqs. ( p9|) and fl32| ) become identical, i.e. the /^-function of 
is continuous (cf. figure p. 

We have argued in a separate paper [[II] that this transition at zero temperature 
may be visible in electron scattering, where real photoproduction with virtuality 
Q 2 = corresponds to a heat kernel cutoff k 2 = and inelastic scattering with higher 
Q 2 to the equivalent k 2 scales. The transition with increasing k from the constituent 
quark to the parton in deep inelastic scattering is a natural consequence in the linear 
sigma model. The vacuum expectation value <pk stabilizes at small values of k and 
the evolution ends with lim^o^fc = fn- When the heavy particles, the sigma 
and quarks, have decoupled, the change of the vacuum expectation value becomes 
proportional to k, the evolution stops (see figure ||). It is interesting to investigate 
how a more sophisticated calculation with running Yukawa coupling gj, and wave 
function renormalization relates the field theoretic model to more observables. 
Thereby the evolution itself may be tested in a region which perturbative QCD 
i. e. the Altarelli Parisi evolution equations cannot reach. 



3.2 The evolution for finite T 

The extension of the flow equations to finite temperature for thermal equilibrium 
is achieved by the Matsubara technique. We integrate the momenta in equation 
(|16|), by splitting the zero momentum component from the three-dimensional spatial 
momentum components and convert the integration over q into a summation over 
Matsubara frequencies uj n for the bosons and v n for the quarks M: 



(0' 



An 2 n 2 T 2 



(2n + 1) 2 tt 2 T 2 



n e 



(33) 
(34) 



In the symmetric phase the corresponding equations are written in terms of the 
positive mass parameter m\ and A*,: 



dvk 
k dk 

k dm 2 
2~dk~ 



k d\ T k 
2 dk 



k 4 7iT 
2(4vr) 2 l r 

k k 2 



(4tt 

ANc9 \ 2 



(47T) 

12(Ap 2 
(4tt) 2 

SN c g 4 
(4tt) 2 



4 E 

. n=— cx 

3vrT 
~2T, 

3ttT 
Y~k 



;i + u 2 /k 2 ) 3 / 2 



{l + u 2 Jk 2 f/ 2 



,(35) 



E 

,=— oo 
oo 

E 



(1 + + ml) /k 2 f 2 
1 



15ttT ~ 
15ttT ~ 



;i + v 2 jk 2 f 2 
i 



(36) 



(l + (ui 2 +m 2 )/k 2 ) 7/2 
1 

(1 + </k 2 ) V2 



(37) 
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The flow equation for the free energy density ffc(O) in the symmetric phase, 
eq. Q55D, does not depend on any coupling neither at zero nor at finite temperature. 
The minimum of the effective potential yields an integration constant. Each contri- 
bution to the flow equations at finite temperature is taken into account in a similar 
way as in the zero temperature case (see above). The modifications in the flow 
equations appear in the threshold functions within the squared brackets in eqs. fl35D 

-(0)- 




L: I I I I I I I I I I 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

T/k 

Figure 3: The ratio of bosonic threshold functions for tempera- 
ture T = and T ^ with different mass parameters y — 0,2,4 
as function of T/k. The dashed line is the function 157r/8 • T/k 
which demonstrates the linear behaviour of the threshold function 
for large T/k in the case of y = 0. 



For the broken phase we get the following expressions: 
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In the limit <frk 0, i.e. at the chiral symmetry breaking scale k x sB the equations 
for the symmetric and the broken phase are continuous. 
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Figure 4: The ratio of fermionic threshold functions for T = 
and with different mass parameters y = 0, 2, 4 as function of 
T/k. 



One sees that the equations arise from the same one loop diagrams as the T = 
equations and therefore depend on the coupling constants in a similar way. Due to 
the three dimensional momentum integrations fractional powers arise in the thresh- 
old functions. 

In the limit of low temperatures T/k —>■ we regain our old expressions at T = 0, 
i.e. we have the relations: 

3ttT ~ 1 1 

^rhoo (i + {;|}/f + ^) 5/2 ~* (i + i/ 2 ) 2 

It is possible to analytically relate all Matsubara sums to the corresponding four 
dimensional integrals in the low temperature limit cf . appendix [A] and ref . ]TB| . The 
above equations guarantee the right matching of the finite temperature equations 
to the zero temperature equations, i.e. in the limit T — > the sets of equations for 
finite and zero temperature become identical. 

We plot in figure |3| the ratios of finite temperature bosonic threshold function 
appearing in the /3-function for over the zero temperature /3-function for 
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Figure ^ shows the equivalent ratios of threshold functions for the fermions. For 
small ratios T/k up to 0.1 the threshold function ratio is constant which means that 
for temperatures small compared to the infrared cutoff scale all Matsubara modes 
are important. The summation is effectively a continuous integration. 

For large ratios T/k the bosonic threshold functions increase linearly in T/k. (cf. 
figure §). This increase is due to the n = Matsubara mode in the frequency sum. In 
the large temperature limit the (3+l)-dimensional system reduces to a 3-dimensional 
system. Thus the phenomenon of dimensional reduction is automatically included 
in the bosonic threshold functions. In the fermionic case there is no zero Matsubara 
frequency, therefore the fermionic threshold functions decrease with T/k and become 
zero. This means that for high T/k the fermions decouple from the further evolution. 
There also exists a stable plateau in the vicinity of the origin for the Fermi -case as 
in the ref. |§. 

The dimensional reduction manifests itself in the critical behaviour, as will be 
shown in the next section. 



3.3 Critical behaviour and exponents 

The idea underlying the finite temperature calculation is the following: At large 
ultraviolet scale k = A and finite temperature T the couplings of the effective theory 
are identical to the T = case, since the finite temperature only modifies the 
boundary in the imaginary time direction but not the dynamics on a small scale in 
space time. This holds as long as the ratio T/k < 0.12, i.e. in the region where the 
ratio of the threshold functions of the ultraviolet theory does not deviate strongly 
from unity (cf. figures §,[§). Thus temperatures smaller than roughly 140 MeV will 
not change the T = couplings at finite temperatures at the ultraviolet scale A = 1.2 
GeV. 

Therefore we can use the same starting parameters as in the T = theory to 
obtain finite temperature results and solve for each fixed temperature T the evolution 
equations as functions of k. The numerical results show the following behaviour: 

With increasing temperature the mass parameter decreases slightly more 
slowly towards the condensation point (cf. figure [5]), but the main effect of the finite 
temperature occurs below k ~ 800 MeV. 

The threshold functions contain an additional damping due to the Matsubara 
frequency. Below the condensation point the boson condensate 4> 2 k reaches less high 
values at finite temperatures than at T = 0. It decreases with k — > for all tem- 
peratures. This effect is due to the pion fluctuations. The critical temperature is 
reached when the condensate tends to zero. The exact value of the critical temper- 
ature depends on the interplay between the relevant sigma mass and quark mass in 
the intermediate k— region. If the "squared sigma mass" 2A&0I is small relative to k 2 
the sigma fluctuations drive symmetry restoration together with the pionic fluctua- 
tions. Our values of the relevant parameter 2\ k (fil at T = are typically around four 
times the quark constituent mass squared evaluated between 100 MeV < k < 300 
MeV. The sigma mass itself is well defined for all k in the case of explicit symmetry 



breaking. The particle data book [^2] gives for the fa- or cx-meson a mass range be- 
tween 400 MeV and 1200 MeV. From the nucleon nucleon interaction, however, we 
know that the intermediate range attraction comes from an equivalent sigma mass 
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Figure 5: The VEV 0^ for k < k x ss{T) and the mass for 
k > k xSB {T) for different temperatures. (Upper line T = 5 MeV, 
middle line T = 150 MeV and bottom line T = 160 MeV). The 
decrease of the order parameter <pk towards zero for k = signals 
the chiral phase transition. 



which has a value of ~ 560 MeV e.g. in the Reid soft core potential p3 |. We add 



as a note of caution that even if the 0(4) dynamics determine critical indices at T c 
correctly, the value of T c may be influenced by gluonic degrees of freedom exterior 
to the linear cr-model in its present form. E. g. one can imagine that the a cou- 
ples to the glueball which dissolves at the critical temperature if the deconfmement 
transition occurs together with the chiral phase transition. With the linear sigma 
model as given and the above parameters we obtain a critical temperature: 

T c « 164 MeV . (42) 

This temperature is higher than the critical temperature of the Wetterich group, 
even if f n = 93 MeV had been used in || for the chiral limit. Wetterich et al. ob- 
tained a critical temperature T c « 100 MeV. They also use different non-analytic 
threshold functions and wave function and coupling constant renormalization. But 
we do not expect the latter modifications to change the critical temperature drasti- 
cally. 

In figure § we show the order parameter <pk=o{T) as a function of temperature. 
At a first glance one might think that the calculated behaviour follows roughly the 
chiral perturbation theory until T « 110 MeV, then it deviates because of strong 
massless mesonic fluctuations which are not properly taken into account in chiral 
perturbation theory. In detail our calculation coincides with chiral perturbation 
theory up to temperatures about 45 MeV then small deviations are visible. In chiral 
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Figure 6: The order parameter <pk=o as a function of the temper- 
ature T. 



perturbation theory the temperature dependence of the light quark condensate for 
massless quarks, plotted in figure |6|, is given by the expression 

mr , T'< I* T> ln A, + 0(T8) (43) 



8/. 2 384/^ 288/| T 
with A q = 470 ±110 MeV Ul4 |. In a model without quark confinement like the one 



used here the fermion loop may be overestimated cf. ||15|| . In general the mesonic 
modes have Boltzmann suppression factors oc exp(— 2uiq/T) compared to a Boltz- 
mann factor cx exp(— rnq/T) for the quark degrees of freedom, where rriQ is the 
constituent quark mass. Therefore the mesons probably contribute to the pressure 
at somewhat higher temperatures than unconfined quarks. 

The renormalization group calculation gives a non mean field behaviour for the 
order parameter at T c , since the long range fluctuations of the softening a modes 
are gradually taken into account by the successive lowering of the cutoff. The renor- 
malization group flow equations at finite temperature sum the diverging thermal 
fluctuations properly near the critical point. 

In the vicinity of T c we obtain scaling behaviour of the order parameter with a 
critical exponent (3 ~ 0.40 as shown in figure In fact it is possible to determine the 
exponent (3 by analyzing numerically the power law behaviour of the order parameter 
near T c . We plot log((f> k ) versus log((T c — T)/T c ) in figure 0. The data points lie on 
a linear curve, thus 4>k scales like 

rp _ rp /3 

A oc (44) 
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with a critical exponent (3 ~ 0.40 which is in good agreement with lattice results for 
the 0(4)-theory in three dimensions [I]. We also compare the result with the mean 
field value 3 = 0.5 in figure [7|. 
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Figure 7: Determination of the critical exponent /3 from the tem- 
perature dependence of the order parameter 0^. 



At finite temperature the limit T/k —>■ oo enforces dimensional reduction. One 
knows, that the linear sigma model lies in the same universality class as the three- 
dimensional 0(4)-Heisenberg model, since the fermions are not contributing to the 
critical fluctuations near T c |13|] . Thus we deal with an effectively three-dimensional 



bosonic 0(4)-symmetric sigma model. For the 0(4)-theory there are six common 
critical exponents but only two of them are independent because they are related 
by four well-known scaling relations of the three-dimensional scalar 0(4)-model: 

a = 2 — dv , 

P = \{d-2 + n), 

7 = (2-77)1/, 
d + 2-r] 



6 



d-2 + r] 



where d is the dimension of the system (d = 3) [p4j] . Using these relations our results 
for 8 = 0.40 and 77 = (no wave function renormalization is used) predict the values 
in tab. ([!]) which are denoted by RG. 

Because we know the free energy density of the system we can calculate the 
specific heat directly , which is presented in the next section and verify the above 
predictions by a determination of the critical exponent a. We find numerically 
a = —0.39 which is consistent with the scaling relation prediction and the 0(4) 
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Mean Field 


RG 


0(4) Lattice 


QCD Lattice 


1/(35 


2/3 


0.5 ±0.06 


0.537(7) 


0.77(14) 


1/5 


1/3 


0.2 


0.2061(9) 


0.21 ...0.26 


1 - 1/5 


2/3 


0.8 


0.7939(9) 


0.79(4) 


(1-/3)/ (35 


1/3 


0.3 ±0.06 


0.331(7) 


0.65(7) 


a//36 





-0.2 ±0.1 


-0.13(3) 


-0.07 ...±0.34 



Table 1: Critical exponents of the three dimensional 0(4) model in mean field, 
in our approach with numerical errors (RG) and lattice calculations compared to 
lattice QCD calculations with staggered quarks (QCD Lattice) (cf. ref. n). 



scaling is confirmed (cf. fig. H and tab. ([!])). One should note though, that the errors 
quoted for our results correspond to numerical errors but do not capture truncation 
errors. 

An especially interesting point is to find the window for the critical dynamics. 
Further investigations away from the critical point may give an indication how far 
away from T c the renormalization group improvement is important for the behaviour 
of the order parameter. In the superconducting phase transition only the mean field 
behaviour is experimentally relevant. Figure ^| shows that in chiral QCD the critical 
behaviour influences a wide range of temperatures, i.e. from T = 130 MeV to the 
critical temperature T c . In real QCD, however, the finite quark masses spoil the 
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second order phase transition. Also gluon effects will change the behaviour of the 
pressure in comparison with the effective linear cr-model. Our findings for the critical 
exponents show excellent agreement with numerical simulation results of the lattice 
0(4)-model. 

Numerical simulations of lattice QCD, however, indicate a critical behaviour of 
the order parameters which is still ambiguous. The critical index 8 gives a strong 
hint for dimensional reduction. For the exponent (3 a possible overlap exists. For 
the specific heat the QCD gluonic degrees of freedom may alter the behaviour sig- 
nificantly. The treatment of light quarks and pions with large correlation lengths is 
problematic on the lattice. 



3.4 Free energy at finite temperature 

Now we turn to the discussion of the equation of state for the linear cr-model in the 
chiral limit. We emphasize this point, because the physics of the model calculation 
is clarified by the equation of state. If the model is complete in the low virtuality 
region, the transmutation of quark-gluonic degrees of freedom occurs above the 
ultraviolet scale. Is this assumption justified? 

Considerable care has been devoted to the treatment of critical fluctuations. In 
order to mark the progress achieved so far, the equation of state is the next goal. Of 
course one does not expect to see free gluons above T c , but what kind of interaction 
effects among the quarks do we have? Can they be compared to deviations of the 
lattice results from the Stefan Boltzmann limit? 

The most important thermodynamical tool for the analysis of a phase transition 
is the free energy, which is given as a logarithm of the partition function. The 
logarithm of the partition function is equal to the effective action, defined in eq. (|T5D . 
Thus the solution of the flow equations, (^) and (|38D , yields the free energy density 
/ as function of the temperature so that the pressure is given by 

P(T) = -v^o = -f. (45) 

The internal energy density is defined by 

u = — (46) 
dT ' K > 

with (3 = 1 /T and the entropy density by 

• - 4 ■ w 

These quantities are connected by the relation P = Ts — u. 

Finally, the specific heat is given by a further temperature derivative of the free 
energy density 

2 du 

= % • (49) 



19 



In order to find the temperature derivative of the free energy density one has 
to take into account the temperature dependence of the quartic coupling Afc(T) and 
the VEV (j>k(T). In principal the thermodynamic functions can be obtained by 
numerically differentiating the k = value of the free energy density as a function 
of temperature. In practice the accuracy of this purely numerical differentiation 
is not sufficient, especially around the critical temperature. We therefore prefer to 
create for each phase a set of additional six evolution equations which contain the 
first and the second derivatives with respect to temperature of the original evolution 
equations, cf. eqs. (ftfi|),(|T7D for the symmetric phase and eqs. (TO) for the broken 
phase. Thus for each temperature derivative one has to solve three additional flow 
equations in this truncation schema, i. e. in total we have nine highly, non-linear 
coupled flow equations. Because they are lengthy and are given by a straightforward 
calculation we renounce to quote the result here. 
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Figure 9: Pressure P/T 4 



Figure 10: Entropy density s/T 3 
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Figure 11: Internal energy density u/T 4 



Figure 12: Specific heat c/T 3 



The procedure to solve these equations is the same as before. The thermody- 
namic predictions are restricted to a temperature window, which is dictated by the 
omission of the temperature dependence of the initial values for the flow equations 
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at the ultraviolet cutoff scale (cf. discussion sec. |3lJ). From this point of view a large 
scale to start the evolution is desirable in order to reach predictive power at least up 
to the critical temperature. For all calculated thermodynamic quantities we choose 
a ultraviolet cutoff A = 1.52 GeV, which is close to the Landau pole of the quartic 
coupling in the chiral limit. Due to this Landau pole inherent in the flow equations 
and the necessity to satisfy certain phenomenological conditions discussed below a 
higher choice for the ultraviolet scale is not possible. We also changed to 88 MeV 
in the chiral limit, so that we obtain for the finite pion mass case the measured pion 
decay constant f n = 93 MeV automatically (cf. sec. ||). With the new ultraviolet 
cutoff we need new initial values. We take Aa = 225, which is larger than the one 
used before with a smaller cutoff due to the growth of the quartic coupling in the 
ultraviolet. We have to adjust the starting mass to m\ = 1 MeV, because for high 
cutoff with large A& the evolution equation of mf (cf. eq.(^)) is governed by the 
bosonic term proportional to A& which first drives the mass up for Ak < before 
the fermion loop decreases again. 

For the flow equations we find for the critical temperature the smaller value 

T c = 149MeV , (50) 

which follows from the smaller but is still much higher than the value found by 
the Wetterich group 0. 

The initial values for the additional flow equations at the ultraviolet scale A are 
zero. At the chiral symmetry breaking scale k x sB ~ 0.8 GeV which is also a function 



of temperature the slopes of the flow equation for in the symmetric phase (p5| ) 
and in the broken phase (|38| ) are identical. The evolution of the potential with 
respect to k already stabilizes at small but finite k— values (k « 150 MeV). 

The result for the normalized pressure (cf. eq. (f£5|)), scaled with the temperature 
P/T 4 , is plotted in figure |9| and in figure |10] the scaled entropy density s/T 3 (crosses) 
is shown. In figure |TI| the scaled internal energy density u/T 4 and in figure [T2| the 
scaled specific heat c/T 3 are shown for the chiral limit. 

In order to compare our results with a gas of non-interacting quarks, antiquarks 
and pions, we evaluate the corresponding thermodynamic functions with the tem- 
perature dependent masses evolved to ~ (dashed lines in the figures). 

For small temperature the quarks are massive and therefore thermo dynamically 
suppressed by their Boltzmann factor. The main contribution to all presented ther- 
modynamic observables in that temperature limit stems from the massless pions. 

We stop the flow equation at a small infinitesimally finite scale k and obtain 
a finite sigma meson mass due to a finite quartic coupling. We have checked that 
there is no influence of a sigma degree of freedom to the thermodynamic functions at 
finite temperature, because the relevant quantity for decoupling is the mass squared 
over k 2 . 

Due to the linear evolution of the quartic coupling with k for k <C T which 
tends to zero for k = in the chiral limit, the sigma meson decouples from the low 
momentum mode thermodynamics and we obtain the expected Stefan-Boltzmann 
limit of a massless pion gas for small temperatures. 

The respective thermodynamic functions for three pionic degrees of freedom are 
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Figure 13: Specific heat c/T 3 in the vicinity of the critical temper- 
ature. 



given by 
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In each figure |9|-|r2| the lower horizontal lines represent the pion gas limit. With 
increasing temperature all observables except the pressure follow roughly the ideal 
gas behaviour. At the critical temperature the quarks are massless and the Stefan- 
Boltzmann limit for a massless quark- ant iquark and massless pion-sigma meson gas 
is depicted by the upper horizontal line in the figures. The thermodynamic functions 
for a gas with N c Nf massless quarks can be calculated analytically: 
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For the numerical calculations at finite temperature one can not simplify the 
computer effort. In ref. the authors reduced the complicated and time intensive 
calculations by substituting the simpler zero temperature flow equations for the 
finite temperature flow equations as long as k > 10T. For smaller scales they used 
the finite temperature flow equations. This is certainly a good approximation for the 
quantities considered there. But the low temperature equation of state would not 
be described correctly in such an approximation. We therefore use during the whole 
evolution the finite temperature flow equations and optimize the convergence of the 
Matsubara series by an generalized 0-function transformation, which is described 
in the appendix. In practice we take more than N = 4000 Matsubara frequencies 
into account in the low temperature region (k ^> T) in order to adequately solve 
the equations over the hole scale range. For higher temperatures (k <C T) only a 
few Matsubara modes are important. 

On the other hand our ultraviolet scale might be identified with the scale 0, 
where the meson bound states appear. Due to the neglect of the wave function 
renormalization we do not suppress the kinetic terms of the mesons at the compos- 
iteness scale. This has an important consequence: Fixing the Yukawa coupling g at 
the compositeness scale to its infrared value which yields a finite constituent quark 
mass of roughly m q « 300 MeV, we do not encounter the partial approximate in- 
frared fixed point, found in 0], which would weaken the dependence of the infrared 
physics on the chosen initial values. In our case we can choose a high ultraviolet 
scale and are nevertheless able to adjust the infrared quantities like e. g. the pion 
decay constant. 

To improve the finite temperature results we use a so-called generalized O- 
function transformation which accelerates the convergence of the Matsubara sums 
leading to modified Bessel function K u (x). This reformulation has the pleasant side 
effect to separate analytically the thermal fluctuations from the T = fluctuations. 
The right hand sides of the flow equations split into the pure zero temperature flow 
equation part and a finite temperature part. We are aware that the separation of 
these terms does not decouple the quantum from the thermal fluctuations because 
of the nonlinearity of the evolution. 

This property, however, allows a clean normalization: neglecting the finite tem- 
perature part of the corresponding flow equation yields the correct zero temperature 
contribution. 

In our approximation of neglecting the scale and temperature dependence of the 
Yukawa coupling as well as the quark and meson wave function renormalization, the 
constituent mass follows the temperature behaviour of the order parameter 0fc(T). 
In the chiral limit above the critical temperature of T c rs 150 MeV the quark masses 
are identically zero and the mesonic masses are degenerate and increase with temper- 
ature because the mesonic coupling also increases with temperature in the symmetric 
phase. For finite pion masses this behaviour is smoother, but qualitatively similar 
and will be discussed in the next section. 
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4 Equation of state for finite pion mass 



It is straightforward to derive thermodynamic quantities for finite pion masses by 
introducing an explicit symmetry breaking term c in the potential. The explicit 
symmetry breaking constant c is kept fixed and is proportional to the current quark 
mass. 

We find the following mass identifications: 

pions: m* >fc (T) = c/</> k (T) (59) 
sigma: m^(T) = c/<p k (T) + 2A fc (T)0 2 fe (T) (60) 
quarks: M q , k (T) = g<p k (T) (61) 

It is not clear if the cr-mass, given above, is actually the physical cr-mass and necessar- 
ily a pole in the two-point Greens function. We follow ref. || where this expression 
is also denoted as the cr-mass and corresponds to the radial mode of the 0(4) model. 
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Figure 14: Mesonic masses versus temperature. 

We fix the constant c via the pion mass (m n (T = 0) = 135 MeV) to c = 1.7 • 10 6 
MeV 3 , which results in a sigma-mass of about 550 MeV at vanishing temperature. 
In reference the sigma mass is considerably smaller. For finite pion masses the 
two sets of flow equations for each phase in the chiral limit merge together to one 
set of flow equations. The constituent quark mass follows the order parameter and 
never becomes zero in the symmetric phase for high temperature due to a smooth 
crossover between both phases. The evolution equations for finite pion masses have 
been studied for two different cutoffs A = 1.2 GeV and A = 2.0 GeV. Due to the 
symmetry breaking term in the potential the initial values are not so sensitive to a 
change of cutoff. The respective initial values are for A = 1.2 GeV A = 115, m = 1 
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MeV and A = 430, m — 1 MeV for A = 2 GeV. In fig. [14| the resulting mesonic 
masses are shown as function of temperature. The scalar meson mass drops relatively 
quickly with temperature while the pseudoscalar masses are almost temperature 
independent in the broken phase. We think that a wave function renormalization 
correction for the scalar fields does not change this situation dramatically at least 
in the broken phase. 
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Figure 15: The order parameter <pk=o for different pion masses. 
(Lower curve: chiral limit; middle curve: m n « 45 MeV; upper 
curve: m n « 135 MeV). 



In general the derivative of these masses with respect to temperature is inde- 
pendent of the choice of the initial UV cutoff. The increase of the masses with 
temperature in the restored phase, however, depends slightly on the choice of the 
UV starting cutoff. We consider this as an artefact of the omission of wave func- 
tion renormalization, which would influence the temperature increase of the quartic 
coupling in the restored phase. 



In fig. ITB we compare the order parameter for different pion masses. The lower 
curve represents the chiral limit where one sees a second order phase transition. 
This transition is smeared out if one switches to finite pion masses. The middle 
curve is achieved for a pion mass of roughly 45 MeV and the upper curve shows 
the realistic case. Note that the order parameter at zero temperature increases with 
the pion mass automatically. If one defines a pseudocritical temperature analog to 
ref. [§] as the inflection point of the order parameter one observes an increase of the 
critical temperature with the pion mass from T c w 150 MeV for the chiral limit up to 
'T' c ~ 190 MeV for a pion mass m n = 135 MeV. Due to the finite symmetry breaking 
term in the potential the order parameter never vanishes at high temperature. 

Above the critical temperature of the chiral transition about 200 MeV we see a 
parity doubling of the scalar and pseudoscalar particles, which is a signal of chiral 



25 



restoration and the masses for the mesons rise proportional to temperature due to 
the increase of the effective mass term with temperature and the symmetry breaking 
term. This behaviour is similar to the behaviour of meson correlation functions for 
massless free quarks above T c on the lattice PB| . 

The hadronic screening masses at high temperature are controlled by the lowest 
Matsubara frequency of quarks and increase therefore like 2ttT. The physics in our 
case is probably more related to the meson dynamics. 
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Figure 16: Scaled pressure (-P(T) — P(T = 0))/T 4 for different 
systems as function of temperature (see text for details). 



In the remaining section we focus on the thermodynamic. In general the follow- 
ing physical scenario becomes evident: For low temperature the quarks and mesons 
are massive, while in the partition function with unbroken chiral symmetry the pions 
remain massless^]. In the chiral limit for low temperature one expects and indeed 
obtains that the physical observables are dominated by a massless non-interacting 
pion gas and that the Stefan-Boltzmann limit is saturated. On the other hand, for 
high temperature in the thermally restored chirally symmetric phase the quarks are 
exactly massless while both, the pions and the sigma mesons, have finite (degen- 
erated) masses due to the restoration of chiral SU{2) x SU{2) symmetry. Thus 
all thermodynamic quantities at high temperature are dominated by the massless 
quarks with N = 4N c Nf degrees of freedom modified by their interaction. 

A similar but not identical situation holds for finite quark masses. The low tem- 
perature gas is hadronic made out of pions mostly, the high temperature plasma is 
a quark plasma. At the critical point the sigma mass is small. Therefore, except 

3 In this case the so-called sigma meson mass, which is proportional to the quartic coupling 
tends to zero in the IR region. 
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near the critical temperature the physical sigma meson contribution to the thermo- 
dynamic functions for low and high temperature can nearly be neglected. In the 
following figures we show the pressure P/T 4 (cf. fig. [16D and the scaled internal 
energy density u/T 4 (fig. |T7p and the scaled entropy density s/T 3 (fig. [PS]). In all 
figures we compare the results with a non-interacting quark-meson gas (uppermost 
curve) and also with a ideal pion gas (lowest curve). 



energy density u/T 4 



ideal quark + meson gas 



A = 2 GeV 




100 150 

Temperature [MeV] 

Figure 17: Scaled energy density u/T 4 for different systems as 
function of temperature (see text for details). 

For finite masses the pressure for small temperature is exponentially suppressed 
as expected. At T ps 200 MeV the total pressure caused by interacting pions, sigma 
and constituent quarks reaches the maximum at (P(T) — P{T = 0))/T 4 = 1.46, 
decreases in the following a bit and settles down around T w 300 MeV at 1.37 (see 

fig- HI). 

A pressure due to the quarks exists also for small temperature, which is a con- 
sequence of the fact that the quarks are not confined in the model. For zero tem- 
perature there is no pressure because the quarks are massive but with increasing 
temperature, the pressure rises and at temperatures around T pa 70 MeV the quark 



pressure is comparable with the mesonic pressure. In fig. [16| the lowermost curve 
shows the pressure of A^ = 3 pions. 

With increasing temperature the quark contribution takes over and the total 
pressure saturates at ~ 55% of the Stefan-Boltzmann limit for T > 250 MeV. This 
result is comparable with an NJL model calculation for small temperature using a 
RPA approximation including fluctuations. Even below the critical temperature a 
strong cutoff dependence enters into the results of previous calculations [^9 . 



At very high temperatures T > 180 MeV our model is no longer adequate: the 
pressure depends on the ultraviolet cutoff which suppresses high momenta. There- 
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Figure 18: Scaled entropy density s/T 3 versus temperature for 
different systems (see text for details). 

fore all non-universal thermodynamical quantities have a tendency to be too small 
at too high temperatures. In order to demonstrate this property we have plotted two 
curves in the figures. The lower line displays the pressure for an ultraviolet cutoff 
A = 1.2 GeV while the upper curve is calculated with A = 2.0 GeV. One sees nicely 
which temperatures can be described for a given cutoff. Both curves coincide up to 
a temperature of T pa 120 MeV, the lower curve turns off for higher temperature. 
Indeed it should be valid only up to T pa 120 MeV. We have checked that the results 
in the shown temperature region are stable with respect to a further increase of the 
ultraviolet cutoff beyond A = 2.0 GeV. 

The uppermost curve in figure [16] displays the total contribution of a non- 
interacting gas of quarks and mesons calculated with temperature dependent masses, 
which are extracted out of the flow equations. It reaches almost the Stefan-Boltzmann 
limit (P/T 4 = 2.7) around T pa 230 MeV due to the decrease of the vacuum ex- 
pectation value and therewith decreasing quark masses. On the other hand, the 
pions and sigma-meson pressure increases im a decreases) up to T pa 200 MeV and 
settles down in the symmetric phase on a constant value. This reflects the fact, that 
the meson masses increase proportional to the temperature and the ratio mass over 
temperature is almost constant. The pressure from the flow equation is considerably 
smaller than this limiting Stefan Boltzmann Plasma. 

In fig. |T7] the energy density, scaled by T 4 , is shown as function of the temper- 
ature. The interpretation of the results and notation of this figure is analogous to 
the pressure calculation. In the chiral limit at low temperatures we observe again 
the Stefan-Boltzmann limit u/T 4 = 7r 2 /10 of the massless pion modes (N = 3). For 
finite pion mass at high temperatures the dependence on the ultraviolet cutoff is 
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slightly stronger compared to the pressure calculation. 

The entropy density, scaled by T 3 , is plotted as function of temperature in fig. |TJ 
and is obtained as a derivative with respect to the temperature of the free energy 
density, s = —df/dT. The low and high temperature features of this quantity can 
be understood in a similar fashion as the pressure and energy density. For finite pion 
masses the entropy density is continuous in the vicinity of the chiral critical temper- 
ature (T c ~ 150 MeV) and reaches 60% of the Stefan-Boltzmann limit around these 
temperatures, which equals to s/T 3 = 8N C ■ 31vr 2 /720 + 16tt 2 /90 = 11.95. Due to 
finite meson masses which increase proportionally to temperature in the symmetric 
regime, the bosonic contribution to all thermodynamic quantities is strong sup- 
pressed. All shown thermodynamic quantities are reliable in a temperature window 
of < T < 180 MeV for an ultraviolet cutoff A of 2.0 GeV. 



5 Conclusions 

We have presented a new method to solve field theoretic problems at finite tempera- 
ture. It is based on a nonperturbative application of the renormalization group. By 
successively integrating out modes in Euclidean momentum shells we can handle the 
nonlinear field couplings in a nonperturbative way, especially when long wavelength 
modes generate singular infrared behaviour near the critical temperature. 

We apply a cutoff function in Schwinger proper time which combines a virtuality- 
or mass- like exponential cutoff with a polynomial which makes the calculation of 
the derivative of the effective potential with respect to the evolution scale insen- 
sitive to the ultraviolet modes. This reflects the modern view of renormalization 
theory, all field theories are effective theories with an ultraviolet cutoff and we al- 
ways follow their infrared behaviour. The evolution equations are finite and well 
behaved, which a simple mass type cutoff would not guarantee. The specifically 
chosen heat kernel cutoff function has the additional advantage that all resulting 
evolution equations can be expressed in a very transparent and analytic form. In 
the one loop approximation propagator like functions for the individual particles 
appear with denominators which contain the respective masses depending on the 
effective potential. The renormalization group allows to improve the one loop ap- 
proximation to sum in an optimal way the infrared singular fluctuations due to 
massless boson modes. 

In this paper we have in addition truncated the effective potential to its lowest 
terms in the boson fields compatible with chiral symmetry. This is a simplification 
which is not necessary and can be improved |TB|,[ff7[]. Also the consideration of 



coupling constant and more importantly wave function renormalization evolution 
equations is planned. The inclusion of vector and axial-vector mesons at finite 
temperature to study their mixing is also possible within this method and work in 
this direction is in progress |19| . 



This new method has been applied to the finite temperature phase transition 
from a constituent quark gas interacting with mesons to a massless quark phase in 
the framework of the linear sigma model. The underlying assumption of a reduction 
of the QCD degrees of freedom to an effective hybrid quark meson field theory at 
an ultraviolet scale of about 1 GeV has been analyzed via the resulting equation 
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of state. Of course, qualitatively the pressure is always dominated by the massless 
modes which are the pions at low temperature and the quarks at high temperature. 
But our calculation has highlighted finer features of the finite temperature behaviour. 

The evolution equations decouple the quarks from further evolution when T/k 
becomes too large. This does not mean, however, that the quarks do not contribute 
to the pressure below the chiral phase transition. The renormalization group tech- 
nique cannot mimic confinement when the field theoretic model does not contain 
this dynamical aspect. The equation of state at high temperatures shows consid- 
erable deviations from a free system of massless quarks. In our interpretation of 
the chiral phase transition long wavelength fermion modes really have zero masses 
at T c . Quarks with moderate or higher virtualities are still massive and only at 
temperatures considerably higher than T c the whole virtuality spectrum of modes is 
massless. 

The real strength of the renormalization group lies near the critical temperature 
where the a modes become massless. Here the linear sigma model with quarks is 
entirely dominated by the bosonic sector. The evolution equations reproduce well 
the universal critical behaviour of the model. We have calculated the critical indices 
/3 and a which are in good agreement with lattice simulations of the 0(4) theory. The 
missing wavefunction renormalization correction is known to be small. The evolution 
equations at finite temperature contain the physics of dimensional reduction: For 
T/k > 0.3 the bosonic loop contributions reduce to their zero Matsubara frequency 
part, i.e. they become equivalent three dimensional loops with effective couplings 
which have acquired temperature factors. The fermions decouple in this region of 
T/k. The differences of some of the critical indices to the mean field behaviour are 
small, i.e. (3 = 0.4 compared to (3 m f = 0.5. This does not mean, however, that the 
region around T c where the critical dynamics is important is small in the equation 
of state. 

The phenomenological relevance of the critical behaviour is more limited because 
of the finite quark mass which breaks chiral symmetry explicitly and smoothes out 
all remarkable features of the 0(4) behaviour. In that respect the chiral dynamics 
is mostly important as a theoretical benchmark. It remains to be seen how far the 
critical 0(4) behaviour is seen in full QCD simulations. 

The evolution equations have the unique property to unify different aspects of the 
hadronic system which conceptually and experimentally are widely separated. The 
zero temperature evolution gives a view of the field theoretic system with different 
resolution. The virtuality cutoff k is intimately related to the resolution of a probe 
with which we may study the system e.g. with photons in deep inelastic scattering. 
At high resolution we see a massless partonic system, at low resolution we encounter 
massive constituent quarks. There are not separate hadronic theories for different 
length scales: a high energy, an intermediate energy and a low energy theory, ideally 
they all evolve smoothly into each other. The challenge of hadronic physics lies in 
the fact that the important degrees of freedom transmute in a strongly coupled 
theory, very much different from a weakly interacting theory. Studying the hadronic 
system at finite temperature and/or finite density allows to emphasize different 
modes by varying the external parameter and thereby showing the transition of 
the dynamics explicitly. The evolution equations unify these different aspects of 
QCD which when studied separately look like rather independent developments in 
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hadronic and nuclear physics. 
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A Matsubara Series: 

In this appendix we summarize the low temperature limit of all in this work used 
Matsubara summations. They can be found by a comparison of the finite tempera- 



ture flow equations with the zero temperature equations. The expression of eq. (41) 
can be generalized by 



oo 



1 X ^Q 2 -4... (2a- 3) 1 1 



X ^ {l + Col + y 2 Y ~^ l-3-5...(2a-2)vr(l + j/ 2 )( 2 «- 1 )/ 2 ^ 
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with x := T/k, u 2 := w 2 /k 2 and y 2 := y 2 /k 2 and a = |, | . . . being a fractional 
power. We set (2a — 3) =: 1 for a — §. 

For the derivatives with respect to temperature of the respective Matsubara sums 
we find for the low temperature limit (/ = 1 describes the first derivative and I = 2 
the second one) 
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If one replaces n — > n + ~ in the last equation one gets the corresponding fermionic 
Matsubara sums with the same low temperature limit. In the low temperature limit 
the difference between the fermionic and bosonic Matsubara sums vanishes. 



B Generalized G-function transformation: 

In order to accelerate the convergence of the Matsubara sums appearing in the 
flow equations it is useful to apply a generalized G-function transformation. This 
transformation is based on the following equations: 

For bosonic Matsubara modes and positive integer p we obtainF] 
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(64) 



A proof can be found in |3C 
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with (2p - 1)!! = 1 ■ 3 . . . ■ (2p - 1) and 
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and for the fermionic case 31 
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(65) 



This transformation can be applied to the finite temperature expression of the 
effective potential eq. (|T6|) and the proper time integration can be performed ana- 
lytically resulting in modified Bessel functions K v (x). 

With the following convenient abbreviations (B = bosons and F = fermions) 



and 
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we obtain the flow equation e.g. for the free energy density: 

W l {x^) + B l {x (7 )-%N c F l {x q )\ 

(67) 

Here one sees the separation of quantum and thermal fluctuations: The right hand 
side of the flow equation splits into the zero temperature flow equation part and a 
finite temperature part. 
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